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We present an anisotropic cosmological model based on a new exact solution of Einstein equations. 
The matter content consists of an anisotropic scalar field minimally coupled to gravity and of two 
isotropic perfect fluids that represent dust matter and radiation. The spacetime is described by a 
spatially homogeneous, Bianchi type III metric with a conformal expansion. The model respects the 
evolution of the scale factor predicted by standard cosmology, as well as the isotropy of the cosmic 
microwave background. Remarkably, the introduction of the scalar field, apart from explaining 
the spacetime anisotropy, gives rise to an energy density that is close to the critical density. As a 
consequence, the model is quasiflat during the entire history of the universe. Using these results, we 
are also able to construct approximate solutions for shear-free cosmological models with rotation. 
We finally carry out a quantitative discussion of the validity of such solutions, showing that our 
approximations are acceptably good if the angular velocity of the universe is within the observational 
bounds derived from rotation of galaxies. 

PACS numbers: 98.80.Hw, 04.20.Jb, 04.40.Nr 
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I. INTRODUCTION 

One of the best established facts in observational cos- 
mology is the isotropy of the cosmic microwave back- 
ground (CMB) This high degree of isotropy explains 
the success of cosmological perturbation theory |^ in re- 
producing the spectrum of anisotropies detected in the 
CMB H]. The measurement of these anisotropies, origi- 
nated from primordial fluctuations, has played a funda- 
mental role in the advent of precision cosmology, allowing 
for the first time the determination of several cosmologi- 
cal parameters and the rejection of a large number of 
cosmological models [|| . 

The isotropy of the CMB, together with the apparent 
homogeneity and isotropy of clustering matter, smeared 
out over scales of the order of 100 Mpc [^|, provide the 
main experimental support for the cosmological princi- 
ple 1^ . The spatial homogeneity and isotropy of the uni- 
verse is incorporated in the standard cosmological model 
by using the Friedmann-Robertson- Walker (FRW) fam- 
ily of metrics to describe the spacetime. It should be 
clear, nevertheless, that the cosmological principle is ac- 
tually a reasonable and fruitful hypothesis, rather than 
a proven fact. In order to clarify this issue, at least from 
a conceptual point of view, we want to show that it is 
possible to construct spatially homogeneous cosmologies 
which are anisotropic but still compatible with the ob- 
served isotropy of the background radiation and the mat- 
ter with strong clustering properties. 

Our starting point will be a spatially homogeneous 
spacetime metric of the Bianchi type III Q subject to 
conformal expansion. As shown in Ref. ^ , the condition 
that the expansion be conformal is crucial for the isotropy 
of the CMB. Part of the matter content will be given by 
a two-component perfect fluid describing the (idealized) 
radiation and dust matter that are present in our uni- 



verse. By its own, this perfect fluid cannot account for 
the anisotropy of the spacetime However, we will 
show that this problem can be solved just by including 
an additional matter source consisting in an anisotropic 
scalar field. In this way, we will be able to construct 
an exact solution of Einstein equations with the required 
properties. 

In addition to their important role in inflationary mod- 
els and Brans-Dicke cosmologies the use of 
scalar fields in cosmology has received renewed attention 
during recent years. The observed relation between lumi- 
nosity distance and redshift for type la supernovae (SNe 
la) has supplied strong evidence in favor of an accelerated 
expansion of the universe Jl2| . In order to explain this 
acceleration and fit the SNe la data, cosmological mod- 
els with a new matter component have been proposed. 
This component, called quintessence, can be modeled by 
a light scalar field with a self-interaction potential and 
with minimal coupling to gravity fl^ . Comparison of 
the theoretical predictions with SNe la and CMB obser- 
vations leads to an estimate for the dark energy density 
of this quintessence field that (roughly speaking) is of the 
order of magnitude of the critical density . 

From this perspective, our approach to construct 
anisotropic solutions can be considered as a new ap- 
plication of scalar fields in cosmology. The anisotropic 
field that we will introduce does not really model a 
quintessence component, because its presence does not 
accelerate the expansion of the universe. Actually, as 
we will see, our scalar field produces neither accelera- 
tion nor deceleration, so that it can rather be regarded 
(apart from the anisotropy) as the limit of a quintessence 
contribution when the acceleration vanishes. Much more 
importantly, it turns out that the dark energy density 
associated with this field is similar to the critical density 
of the model at late times. In this sense, our model sug- 
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gests that the inclusion of an anisotropic scalar field may 
provide a mechanism to generate quasiflat |ll5|| universes. 

In the case of our anisotropic cosmological solution, 
the quasiflatness can be rooted in the fact that the en- 
ergy density of the anisotropic field is proportional to the 
inverse square of the scale factor, which is precisely the 
type of dependence that one would expect for the critical 
density at the final stages of the expansion. The reason 
is that, in a Bianchi type III universe, there exists a neg- 
ative contribution of curvature to the energy density, just 
like in an open FRW model, so that at large cosmological 
times one would expect the scale factor to be inversely 
proportional to the Hubble constant and, therefore, to 
the square root of the critical density. 

As we have said, the evolution of the anisotropic solu- 
tions that we will present is not accelerated, so that our 
cosmological model cannot be considered realistic in this 
respect. One could remedy this situation by including an 
additional matter source, given by a quintessence field or 
a cosmological constant. Instead of proceeding in that 
way, we have preferred to keep the model as simple as 
possible, both to isolate the cosmological consequences 
of the anisotropic scalar field and to obtain a solution of 
Einstein equations in which one can complete all calcula- 
tions of cosmological parameters using exact expressions. 

Apart from clarifying the role that anisotropic scalar 
fields may play in cosmology, we will also discuss in de- 
tail another possible application of our exact anisotropic 
solution. This application follows from the fact that 
the spacetime metric of our model is just the vanishing- 
rotation member of a family of spatially homogeneous, 
rotating and shear-free metrics that were studied by Ko- 
rotkii and Obukhov (KO) This family is an ex- 

panding version of a class of Godel-like stationary metrics 
analyzed by Rebougas and Tiomno (RT) [0 . From now 
on, we will refer to them as the RTKO metrics. In the 
limit of small rotation, it is not difficult to employ our 
exact anisotropic solution to construct approximate so- 
lutions for the rotating RTKO cosmologies and supply 
them with a physically acceptable matter content. 

We will also carry out a quantitative analysis of the 
validity of our approximations when a small rotation 
is present. In particular, we will find a bound on 
the angular velocity of the universe guaranteeing that 
our approximate solutions cannot be distinguished from 
the unknown exact ones, at least as far as the energy- 
momentum tensor is concerned. 

The RTKO metrics share in fact most of the good 
properties presented by our exact anisotropic solution 
with vanishing rotation. For instance, they are all of 
the Bianchi type III and possess a conformal Killing vec- 
tor field (CKVF). The existence of a vector of this kind 
parallel to the four- velocity of dust matter and radiation 
is known to be the necessary and sufficient condition for 
the absence of parallax effects and turns out to guar- 
antee the isotropy of the CMB jgf . In addition, the space- 
time does not contain closed timelike curves (CTC's) un- 
less the rotation is considerably large |^. On the other 



hand, it has been shown that the RTKO metrics repro- 
duce the open FRW metric in the limit of small rotation 
and nearby distances All these properties clearly 

make of the RTKO metrics natural candidates to describe 
anisotropic cosmological scenarios with rotation. 

The rest of the paper is organized as follows. Since the 
spacetime of the exact cosmological solution that we will 
construct belongs to the RTKO family, we will first ana- 
lyze the most relevant properties of these geometries and 
obtain their Einstein equations in Sec. II. In Sec. Ill we 
present our exact anisotropic solution with vanishing ro- 
tation, and discuss the corresponding cosmology in Sec. 
IV. Based on the case of zero angular velocity, we find in 
Sec. V approximate RTKO solutions with rotation. We 
also study their cosmological parameters and show that 
the quasiflatness of the cosmological model persists in the 
presence of rotation. By comparing the Einstein tensor 
of our approximate solutions with the energy-momentum 
tensor of its assumed matter content, we derive in Sec. 
VI a bound on the angular velocity of the universe en- 
suring that the relative errors committed in the Einstein 
equations are small. The conclusions of our work are in- 
cluded in Sec. VII. Finally, an appendix containing some 
calculations is added. 



II. THE RTKO METRICS 

The RTKO metrics are described by the line element 

ds^^a'^{ri)[-{dri + le''dyf+dx^ + e^''dy^ + dz^] , (2.1) 

where rj is the conformal time and x, y, z are the spatial 
coordinates, all of them assumed to run over the real axis. 
The parameter /, on the other hand, is a constant that 
can be restricted to be non-negative without loss of gen- 
erality. From now on, we call it the rotation parameter. 
The above spacetime does not contain CTC's if and only 
if I belongs to the interval [0, 1), since it is only then that 
the metric induced on the sections of constant time is 
positive definite In the following, we will restrict 

our considerations to this causal sector, < ^ < 1. 

We will employ the following notation. Greek letters 
will denote spacetime indices, and the indices {0, 1, 2, 3} 
will designate, respectively, the coordinates {?], a;, y, z}. 
In addition, we adopt units such that SttG — c — 1, G 
being Newton constant. 

Metric (|]^) is a spatially homogeneous, Bianchi type 
III metric, with three Killing vector fields given by 

C(i) =d:c- ydy, $(2) = dy, C(3) = dz- (2.2) 

The metric possesses also a CKVF, namely, = 5^. 
Korotkii and Obukhov have proved that, assuming co- 
moving radiation and a comoving observer, the existence 
of a CKVF guarantees the isotropy of the detected CMB, 
with the radiation temperature falling with the inverse of 
the scale factor, and ensures that the redshift of the light 
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coming from astrophysical objects does not depend (ex- 
plicitly) on the spatial positions of the source and the 
receiver, but only on the emission and observation times 
I^Jl^. Besides, the CKVF prevents the appearance of 
parallax effects Furthermore, in these spacetimes 

the shear tensor vanishes for comoving observers (with 
four- velocity equal to u'^ — (5q /a), whereas their rotation 
tensor 0;^,^ is different from zero [0,^. In particular, their 
angular velocity is a; = ^ ujij^^uji^^ jl = I /{2a). 
With the help of the coordinate transformation 



— cosh r + cos sinh r, 
ye^ — sin cj) sinh r, 

rj = fi — l(f> + 21 arctan [ e^"" tan ■ 



(2.3) 
(2.4) 

(2.5) 



one can write metric (2.1) in the cylindrical form |17| 

ds^ = -a^{r]) dfi + 2lsmh^(^-^d(f>^ 

W{r]) [dr^ + sinhV dcj)^ + dz^] . (2.6) 

Here, 77 is the new real time, r is non-negative, and (p is 
an angular coordinate. Note that the scale factor is not 
constant on the sections of constant time fj unless the ro- 
tation parameter vanishes [or a{rj) is a constant number], 
because a depends on the radial and angular coordinates 
r and </> when / ^ . On t he o ther hand, using the change 
of coordinates (2.3) and ( ^.4[ ), one can easily check that 
the sections of constant time rj are the direct product of 
a real line and a two-dimensional pseudosphere. 

Let us finally consider the associated Einstein equa- 
tions. For the diagonal components of the Einstein and 
energy-momentum tensors, one obtains 



riV=(i-^')a^-3(i-;V 



(2.7) 



Tla^ ^ T^a^ = -a^ + (1 - P)a^ ~ 2(1 - P)ad, (2.8 



(2.9) 



where the overdot denotes the derivative with respect to 
the conformal time rj. For the non-diagonal components, 
on the other hand, the result is 



rrl 4 ;2 • 

Jq a = i aa. 



T^a^ = 2/e-^(aa - 2d^), 

T°a^ = (1 - l^)le^a\ 
T^a^ = ~{1 - l^)le^aa. 



(2.10) 



The remaining components of the energy-momentum ten- 
sor must be identically zero. As far as we know, no phys- 
ically admissible matter source has been proposed up to 
date leading to a solution of the above Einstein equations 
when the scale factor is not constant. Hence, no explicit 



RTKO cosmological model has been constructed so far. 
In the next section, we will present an exact solution 
for the case of vanishing rotation that has an acceptable 
energy-momentum tensor. This solution represents an 
anisotropic universe in continuous expansion. 



III. THE ANISOTROPIC SOLUTION 

We will now restrict our attention to the RTKO metric 
obtained when the rotation parameter I vanishes. In the 
absence of rotation, the Einstein equations require the 
energy-momentum tensor to be diagonal. The diagonal 
components must satisfy 
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ea = Aa — a , 
Pio^ — p2a'^ = d^ — 2ad, 



Pia 



1. 



(3.1) 
(3.2) 
(3.3) 



Here, we have adopted the notation e = — Tq for the en- 
ergy density and pi = (z = 1, 2, or 3) for the principal 
pressures of the system. 

Let us start by assuming that ea^ vanishes in the limit 
of infinite scale factor, as would happen if the matter 
content co nsist ed exclusively of dust and radiation |6|j2C 



From Eq. (3J^) one then easily sees that c? ja^ must tend 



to 1/3 when a cxd. Note also that this equation ensures 
that a increases unboundedly with the conformal time, 
provided that d is initially positive. In addition, suppos- 
ing that ea^ is a smooth function of the scale factor, the 
condition ea^ implies that {de/da)a^ vanishes when a 
becomes infinitely large. Using these facts and taking the 
time derivative of Eq. (p.l|), it follows that the quotient 
d/a must also tend to 1/3 when a approaches infinity. 
Employing now Eq. ( |3.2D , one concludes that pia^ has 
a finite limit when a 00. As a result, the dominant 
energy condition ||2|] is violated during the evolution, 
because for sufficiently large scale factors the pressure pi 
becomes larger than the energy density. 

Therefore, if we want to reach an acceptable solution 
of the Einstein equations that respects the energy con- 
ditions, we must include matter sources whose energy 
density does not fall faster than 1 /a^ when the scale fac- 
tor expands to infinity. Probably, the simplest way to 
do this is by introducing an anisotropic massless scalar 
field minimally coupled to gravity. As we will see below, 
the corresponding energy density satisfies precisely the 
minimal requirement of being proportional to the inverse 
square of the scale factor. Furthermore, the inclusion of 
such a scalar field will actually suffice to explain all the 
anisotropics of the model, allowing the rest of the matter 
content to be isotropic. 

In curved spacetime, a massless minimally coupled 
scalar field satisfies the equation 



0, 



(3.4) 
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whereas its energy-momentum tensor has the form 
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(3.5) 



Here, g and g^^ are the determinant and the inverse of 
the four-metric, the semicolon denotes covariant deriva- 
tive, and (5!^ is the Kronecker deha. 

Let us then consider an anisotropic scalar field given by 
$ = Cz, with C being a constant. This kind of source for 
the RTKO metrics was already suggested by Rebougas 
and Tiomno in a stationary c ont ext with rotation p^ . 
It is easily checked that Eq. (^) is in fact satisfied by 
our fi eld in any of the spacetimes (2.1). Besides, from 
Eq. ( |3.5|) , the solution $ = Cz has a diagonal energy- 
momentum tensor, with the following energy density and 
principal pressures: 



'P2 



Here, the superindex (s) refers to the contribution of the 
scalar field. 

Note that the corresponding energy density falls with 
a^, as we had anticipated. In this respect, it is inter- 
esting to note that such a kind of decay for the energy 
density is also expected on the basis of quantum cosmo- 
logy arguments and might even provide a way to solve 
the cosmological constant problem p^ . 

In addition, the principal pressures are now aniso- 
tropic. Using this property, it is actually very simple 
to remove any trace of anisotropy from the Einstein 
equations of our spacetime. Defining e = e -\- e^'^^ and 



Pi= Pi+Pl ', we see from Eq. ( |3.3D that the anisotropic 



contributions of the model can be absorbed in the scalar 
field just by imposing that = 1. Since the orientation 
of z can be inverted at will (producing an apparent flip 
of sign in the cons t ant C), we will fix from now on $ = z. 
Equations ( |3.l| )-( |3^ ) become then 



ea 



pa 



3a^ 



2ad+ y. 



(3.7) 
(3.8) 



where p = pi for any i = 1, 2 or 3. Remarkably, these are 
exactly the Einstein equations of an open FRW model 
with curvature parameter k equal to —1/2 Equiva- 
lently, they can be written as the equations of the stan- 
dard FRW model with k = —1 under the scaling: 



V2r]F 



(3.9) 



Here, the subindex F denotes the conformal time and 
scale factor of the open FRW cosmology. Notice that 
these relations imply that the cosmological time of our 
model coincides with that of the standard FRW space- 
time, because adij = apdr/p. 



From the above comments, it should be clear that the 
evolution of the scale factor in our model reproduces the 
expansion found in an open FRW cosmology, except for 
some qualitatively irrelevant scalings by factors of the or- 
der of the unity. Owing to this fact, and leaving aside the 
anisotropy of the sections of constant time, the cosmolog- 
ical solution that we will construct leads essentially to the 
same history of the universe as a standard open FRW sce- 
nario, at least during the epoch in which the scalar field 
has a negligible contribution to the energy and pressure 
of the system. Furthermore, regarding the anisotropy of 
the spatial sections, we recall that our spacetime metric 
can be written in the form (2^) with 1 = 0. In fact, since 
such metric reduces to an open FRW metric in the limit 
of nearby distances r <^ 1 |1^ , no differences should be 
expected in physical processes or observations which do 
not involve distant regions. 

We are now in an adequate position to obtain the so- 
lution of the Einstein equations that we were seeking. 
In addition to the anisotropic scalar field, we suppose 
that the matter content is given by radiation and dust, 
as it is usually done in standard FRW cosmology. We 
will describe these matter sources by a two-component 
perfect fluid, with comoving four- velocity = ^o/*^- 
The assumption that the radiation present in the sys- 
tem adopts the form of a comoving perfect fluid, together 
with the properties of the RTKO metrics Q , guarantees 
that the CMB of the model is isotropic. Similarly, the 
fact that the dust matter can be treated as a comoving 
perfect fluid ensures the applicability of Hubble law (in 
the leading-order approximation) to any kind of radia- 
tion that could be emitted by dust particles, since the 
radiation frequency varies then just like the inverse of 
the scale factor 1^. For such a matter content, the ex- 
pression of the energy and pressure that appear in Eqs. 
(13) and (O) are 



D 

7^' 



P ■ 



3^' 



(3.10) 



where A and D are two non-negative constants. The 
first term on the right-hand side of these equations cor- 
responds to the radiation component, whereas the dust 
matter contributes only to the energy density |^pO[ | . 



With the above energy and pres sure, Eq. (3.7) turns 



out to be a first integral of Eq. (3.8), and admits a unique 
increasing solution that vanishes at = 0. The exact 
solution is given explicitly by 



D 
3" 



».h,^) 



1 



H^.i.h(^). (3.11) 



This expression can be inverted in 77 > 0, obtaining 
77 = \/21n 



3a + D + \/9a2 + 6Da + 6^2 



d + Vga 



(3.12) 



On the other hand, integrating dt = adrj, we arrive at 
the following expression for the cosmological time: 
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D 



-1] 



2A 

7! 



cosh ( I 
V2J 



(3.13) 



From the last two formulas, one can also calculate t as a 
function of the scale factor. 



IV. THE ANISOTROPIC COSMOLOGICAL 
MODEL 

In the preceding section, we have constructed an exact 
solution of Einstein equations that describes an expand- 
ing universe containing an anisotropic massless scalar 
field and a comoving perfect fluid composed of radia- 
tion and dust matter. The spacetime metric is given by 
the element of the family (2.1) with vanishing rotation. 
As a consequence of the properties of the RTKO met- 
rics, the CMB of the model is isotropic and the redshift 
of the radiation emitted by the comoving dust depends 
only on the emission and observation times |^ . We have 
also seen that our anisotropic metric coincides with the 
metric of an open FRW universe in the limit of nearby 
distances. Moreover, the conformal expansion of our so- 
lution reproduces (apart from some trivial scalings) the 
evolution encountered in a standard open FRW cosmo- 
logy with matter content formed exclusively by isotropic 
dust and radiation. Hence, the history of the universe 
in our anisotropic model parallels that of an open FRW 
solution, at least as far as the scalar field does not sup- 
ply the dominant contribution to the energy-momentum 
tensor. 

Like in the analogue FRW cosmology with energy den- 
sity and pressure given by e and p, the radiation dom- 
inated era of our anisotropic model corresponds to the 
epoch with < a < A^/D. At small times, the universe 
expands from an initial singularity following exactly the 
same evolution law as in standard FRW cosmology (|J2^ , 
namely, 
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(4.1) 



This behavior can be easily obtained from Eq. ( 3.11 ) 
in the region 77 ^ 1. As a particular consequence, 
the Hubble parameter and the energy density adopt, 
at the initial stages of the expansion, the expressions 
H = aja^ = 1/(2*) and e = 3/(4^^), which coincide with 
the result of the standard model in the radiation era. 
In particular, it follows that the initial relative energy 
density is = e/(3i/2) = 1. 

When a increases beyond /D, the dust component 
starts to supply the major contribution to the energy den- 
sity and the universe enters a dust dominated era with an 
evolution of the scale factor similar to that presented in 
an open FRW cosmology. Such era ends when the energy 



of the anisotropic scalar field becomes the most impor- 
tant matter component. This occurs when e*^"^ equals 
the dust energy density, i.e., when a = 2D. We assume 
that A <^ V^D, so that there exists a sufficiently large 
epoch A^/D < a < 2D dominated by matter with strong 
clustering properties. As far as a < 2D, the contribution 
of the anisotropic scalar field is subdominant, and the 
model leads essentially to the same cosmological predic- 
tions as an open FRW model. 

For scale factors larger than 2D, the anisotropic scalar 
field dominates the evolution. The expansion is then of 
the approximate form 



(4.2) 



as one can check from Eq. ( 3.11 ) by analyzing the sec- 
tor of large times. Note that this evolution is linear in 
the cosmological time, like at the final stages of an open 
FRW model. This was in fact expected, because the time 
dependence of the scale factor must always be similar to 
that of an open FRW universe without scalar field, as we 
showed in the preceding section. In the limit rj 00, the 
Hubble parameter displays then the behavior H — l/t, 
and the energy density is e^*-* = H^. Hence, at large 
times, the relative ener gy d ensity becomes f2 = 1/3. 

Actually, from Eq. (3.11) we can derive the exact ex- 
pressions of the Hubble parameter, the deceleration pa- 
rameter q, and the relative energy density Q at all times 
of the evolution. We get 



3a2 + 2Da + 2A'^ 



6a4 
2Da + 2A^ 



3i72 



^ aa 
^ a? 3a2 



3a2 + 2Da + 2A'^ ' 
Da + 2A^ 
+ 2Da + 2^2 



>0, 



(4.3) 
(4.4) 
(4.5) 



where we have used that e = e^'*-' + e and employed Eqs. 



(3.(:) and (3.10). We recall that the parameter q is posi- 
tive when the expansion decelerates. 

In the limits a — s- and a — > 00 (i.e., when rj tends 
to zero and infinity, respectively), we recover from these 
equations the behavior discussed above for H and O. Fur- 
thermore, it is not difficult to prove that Eq. (4.4) defines 
a strictly decreasing function of the scale factor, il{a). 
Since the universe is always expanding in our solution, 
we conclude that the relative energy density of our model 
suffers a continuous decrease from its initial unit value at 
the big-bang singularity, reaching the asymptotic lower 
bound of one-third in the limit of large times. In this 
way, the contribution of the anisotropic scalar field guar- 
antees that the energy density of the model is of the order 
of the critical one during the whole evolution, leading to 
a quasiflat universe. 

Obviously, the model is not fully realistic; in particu- 
lar, the positivity of Eq. (4^) means that the expansion 
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decelerates in our solution, contradicting the present ob- 
servations of SNe la ||l^ . The result q > can be easily 
understood on the basis of our matter content: as we 
have seen, the scalar field leads to a uniform expansion, 
linear in the cosmological time, whereas the presence of 
radiation and dust decelerates the expansion. Note, how- 
ever, that the deceleration is similar to that found in 
a standard open FRW cosmology without cosmological 
constant and quintessence fields. This follows from the 
fact that the deceleration parameter q reflects only the 
time dependence of the scale factor, and this dependence 
coincides in our solution and in an open FRW model. 

In order to attain an accelerated expansion in our 
anisotropic scenario, we could simply add a positive cos- 
mological constant A to the matter content. Indeed, it is 
easy to check that, for t he e poch in which A dominates 
the energy density, Eq. would lead to an exponen- 

tial expansion. Like in standard FRW cosmology, how- 
ever, we have preferred to analyze here the case without 
cosmological constant (or quintessence) because in this 
way we can obtain an explicit solution that allows us to 
perform all calculations to conclusion. In addition, the 
inclusion of other matter sources would have prevented us 
from clearly isolating the consequences of the anisotropic 
scalar field. 

In order to estimate the values of the parameters A 
and D and the present values of a, t, q, and fl in our 
model, we can proceed as follows. From Eq. (O) we get 



a-o 



ml 



2er — 2ed 



(4.6) 



where the subindex means evaluation at the present 
time, and the subindices r and d denote the radiation 
and dust components of the matter content. In addition, 
if f2d is the contribution of dust matter to the relative 
energy density and Z^q is the redshift corresponding to 
the equilibrium between dust and radiation, we have that 
ed = SfidiJ^ and Ir = ed{l + Zeq)^^ . Finally, = e^o^ 



and D — £dO,Q- With these values and formulas (4.4) 



and (4.5), we can also determine the quantities qo and 
flQ. Using (approximately) the values of the concordance 
model 23 1 for the present Hubble parameter and relative 
energy density of pressureless matter, Hq = 65 km/(sMp) 
and fid = 0.35, as well as Z^q + 1 = 5000, we obtain that 
A=1.6x 10^4 m, D = 1.0 x lO^^ m, oq = 1.2 x lO^^ m, 
^0 = 12 Gyr, qo = 0.18, and = 0.57. 

From these estimates, we see that the assumption 
A <C V^D is actually satisfied in our solution. The 
dust era corresponds to the interval of scale factors 
2.5 X lO^'^ m < a < 2.0 X 10^^ m, which is large enough 
for structure formation and contains the present period 
of the evolution. We also see that the equilibrium be- 
tween dust matter and the anisotropic scalar field would 
be reached when a — 2D — 2.0 x 10^^ m, a value of the 
scale factor that is only slightly larger than the present 
one. Thus, in our cosmological model, we would be al- 
most at the end of the dust dominated epoch. 



It is worth noting that, although the CMB of the model 
is isotropic and the redshift of the radiation emitted by 
dust particles depends only on the value of the scale 
factor at the moment of emission, and not on the spa- 
tial position of the source, the fact that the metric is 
anisotropic implies that the distance to astrophysical ob- 
jects with identical redshift varies with the direction of 
observation. One might then worry about the compat- 
ibility of this anisotropy with the available data about 
extra-galactic sources at high redshift, e.g. with the ap- 
parent isotropy detected in the Hubble diagram for SNe 
la at redshifts of order unity. In order to discuss this 
issue, let us consider the angular diameter distance |6], 
which can be defined by the relation dA^ — dflo p4|. 
Here, dAe is the (infinitesimal) intrinsic perpendicular 
area of the source, which subtends the solid angle dilo at 
the origin where, using the homogeneity of the spacetime, 
we locate the receiver ||2^,^. The luminosity distance 
is then r; = ra{l + Z) , with Z being the redshift of 
the source Hence, one only has to care about the 

anisotropics that appear in Tq. Using the expressions 
given in Ref. |2^] (or just applying the formulas of Ref. 
26|), it is possible to show that 



rl = a^iVe) ivo - Vef F[sine'(r/o - ?/e)], 



(4.7) 



where Y{u) = svahu/u, 9 e [0,7r] is the angle formed 
by the line of sight and the z axis, and rj^ is the con- 
formal time of emission. As anticipated, depends 
on the direction of observation and, for fixed Z (and 
present time 770), its maximum tm and minimum are 
reached when sin 9 equals the unity or tends to zero, re- 
spectively. The magnitude of the relative variation of ra 
on the celestial sphere can be described wit h th e quan- 
tity £a = {tm - rrn)lrm. Employing Eqs. (|3.12| ), ^), 
and 1 + Z = aola(ric,), it is straightforward to see that 
£a increases with Z . More importantly, substituting the 
values of the constants A, D, and oq obtained above, one 
can check that the relative variation of the angular dia- 
meter distance is only of the order of 5% for Z = 1, while 
for Z — 2 Ea close to 10%. These variations do not 
seem to conflict with the observational data, and do not 
dominate over the systematic and statistical uncertain- 
ties, evolution effects, and experimental errors that are 
present in the determination of astronomical distances. 

Finally, let us point out that the age of the universe 
in our model (to = 12 Gyr), although very close, is still 
beyond the lower bounds obtained from radioactive dat- 
ing of stars [|7| or studies of globular clusters . These 
results show that (except for the absence of acceleration 
and the corresponding quintessence contribution to the 
relative energy density) our anisotropic model is at least 
compatible with the main features of modern standard 
cosmology. 
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V. APPROXIMATE ROTATING SOLUTIONS 

In this secti on, we will present a generalization of the 
solution (3.11) for a non-vanishing rotation parameter, 
I 0. We will assume the same matter content as in 
the absence of rotation, namely, a two-component perfect 
fluid, formed by radiation and dust, and an anisotropic 
scalar field $ = z minimally coupled to gravity. For small 
values of the parameter we will see that the RTKO met- 
ric that we will obtain can be regarded as an approximate 
solution of the Einstein equations. In this way, one can 
construct an approximate cosmological model describing 
the expansion of a rotating anisotropic universe which 
contains isotropic background radiation. Actually, sup- 
posing that / is sufficiently small, the inclusion of rota- 
tion produces only small corrections in the cosmological 
model constructed in Sec. IV. As a consequence, our ap- 
proximate solutions will lead to a similar cosmology, both 
qualitatively (apart from the existence of an angular ve- 
locity) and quantitatively. 

The energy-momentum tensor will have the form 



(p -t- e)u''u^ 



(5.1) 



where u^^ is the four-velocity of the two-component fluid, 
its energy density and pressure are given in Eq. (3.1C), 
and T^*' denotes the energy-momentum tensor of the 
anisotropic scalar field. Th e components of this diagonal 
tensor appear in Eq. ( |3.6| ) (with C = 1). The parame- 
ters A and D, which determine the energy density, are 
assumed to be exactly the same as in the solution with 
vanishing rotation. Like in that case, we also consider 
comoving perfect fluids with — Sq/u. 

Using the general RTKO non-diagonal metric ( |2.1| ) , we 
obtain the covariant four-v eloc ity = — a(^^ -I- le^Sf^). 
Then, from our deflnition (5J), we see that the diago- 
nal components of the energy-momentum tensor are (for- 
mally) the same as in our solution with zero angular ve- 
locity, whereas all the non-diagonal components vanish 
except T2 . This last component takes the expression 



rjiO 

^2 



4^2 
3^ 



D 



(5.2) 



Let us first consider the diagonal time component (2.7) 
of the Einstein equations. When the rotation parameter 
does not vanish, this equation has the following solution 
for our value of the energy density: 




where we have introduced the definitions 



X. 



1 



(5.3) 



(5.4) 



The above scale factor increases with the conformal time 
in 7/ > an d vanishes at ry = 0. In addition, it reproduces 
Eq. (3.11) when I vanishes. Note also that, since we 
have imposed that / e [0, 1), the ranges of Xi and Y/ are, 
respectively, (1/2, 1] and (0, 1]. 

Substituting the above time dependence of the scale 
factor and t he ex pressi on o f the energy-momentum ten- 
and (E^), it is easy to check that the 



sor in Eqs. (2 



Einstein equation G\ — T3 is satisfied exactly; however, 
the other diagonal spatial components of the Einstein 
tensor differ by a term P /{2a'^) from their assumed va- 
lues. In other words. 



(5.5) 



Concerning the non-diagonal components (2.1C) of the 
Eins tein equations, it is not difficult to prove using Eq. 
( pT^ ) that, when I is small, the Einstein tensor of the 
analyzed RTKO metric provides an approximate solu- 
tion up to terms of the order of P for Gq — Tq and 
of order I for the rest of equations. Therefore, we con- 
clude that the difference between the components of the 
energy-momentum tensor of our system and those of the 
Einstein tensor of the metric (2.1) and (5.3) vanish at 
least as fast as I when I — > 0, and become, in general, 
negligible when the rotation parameter is small. In Sec. 
VI we will use this fact to set an upper bound to the 
global angular velocity in order to ensure that the rel- 
ative error committed in the energy-momentum tensor 
with our approximation is smaller than a certain quan- 
tity. 

Let us now analyze the behavior of our approximate 
cosm ological solutions with rotation. Inverting relation 
( [S.3| ), we obtain the conformal time 




3Xia + D+ ^9Xfa^ + 6XiDa + 6XiA^ 



(5.6) 

and, integrating dt = adr], we get the following expression 
for the cosmological time: 




(5.7) 

From relation (5.3), one can also derive the Hubble 
parameter, the deceleration parameter, and the relative 
energy density of our approximate solutions: 



H 



n = 



2Da + 2^2 



2Da + 2A2) 



3Xia^ + 2Da + 2A^ ' 

Da + 2^2 
iXia? + 2Da + 2A^ 



> 0. 



(5.8) 

(5.9) 
(5.10) 
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These formulas replace Eqs. (|4.3|), ( [4.4D , and (|4.5|), re- 
spectively, when the rotation differs from zero. 

In the limit a ^ 0, we get again H = l/{2t) and q — > 1, 
as in the standard cosmological model. In this limit, the 
relative energy density takes the value Q. = Yi = 1 — P, 
so that e — 3(1 — P)/4t^. The expansion and history 
of the primordial universe is therefore aff ecte d only by 
corrections of the order of P [see also Eq. (5.3)]. On the 
other hand, in the sector of large scale factors a — > oo, 
one can easily check that H = 1/t and g — > 0, just like 
on the exact solution presented in Sec. III. At this final 
stage of the expansion, the relative energy density tends 
to Yi/{3Xi), a limit which is positive for I G [0,1) and 
differs from the value of 1/3, corresponding to the non- 
rotating case, by terms of the order of P, supposing that 
the rotation parameter is small. 

Finally, it is not difficult to prove that the relative en- 
ergy density (5.9) is a strictly decreasing function of the 
scale factor. Like in the model discussed in Sec. IV, 
remains then bounded away from zero during the whole 
evolution, the lower bound being its positive limit when 
a ^ oo. Actually, if / <C 1, the energy density is always 
of the same order of magnitude as the critical one. There- 
fore, we see that the introduction of an anisotropic scalar 
field leads to a quasiflat universe even in the presence of 
rotation. 



sor (namely, T2) is also smaller than the energy density. 
Therefore, with our definition, we are just comparing the 
errors made in the estimation of the energy-momentum 
tensor with its dominant component. 

To analyze these errors, we need to deal with factors 
of the form e^^ that appear in most of the non-diagonal 
comp onents of the Einstein tensor, as can be seen in Eqs. 
( 2.1c ). In doing this, we will proceed as follows. Since the 
model is spatially homogeneous, we can always locate the 
observer at the origin. From a physical point of view, the 
only phenomena that can affect the observer at a generic, 
present time rjo are those that occurred in the spacetime 
region that is causally connected with him. Thus, from 
now on we will restrict our discussion to that region. Let 
us also suppose that we are only interested in events that 
happened in a certain interval of time rj € [rji, tjq] , with 
< ?7i < ?7o- Although we will make rji — at the end of 
our calculations, we prefer to leave this number free for 
the moment in order to allow for other possibilities. 

In a RTKO spacetime, one can check that the maxi- 
mum absolute value that the coordinate x can take at 
time 77 < 770 in the region that is causally connected 
with the origin at present is (770 — t?) / Vl — l^- A point 
at time 77 with this value of x is connected with the 
origin at rjQ by the null geodesic with vanishing z and 
dy/drj = le~^ /{I — P). Hence, the region of the space- 
time that we want to analyze is contained in 



VI. VALIDITY OF THE APPROXIMATION 

In this section we want to carry out a quantitative 
analysis of the error committed in Einstein equations by 
identifying the energy-momentum tensor ( p5.1[ ) with the 
Einstein tensor of the RTKO metric whose scale factor 
is the time function ( |5.3| ). More specifically, we want to 
show that it is possible to set an upper bound to the 
rotation parameter (and hence to the present angular 
velocity) so that the relative error in our estimation of the 
energy-momentum tensor is smaller than a fixed quantity. 

For each component of the Einstein equations, we de- 
fine the relative error introduced with our approximation 
as the quotient |GJ^ — T^|/e, where e = — Tq is the energy 
density of the matter content. We want to analyze under 
which circumstances these relative errors are smaller than 
a given number A. Since, for any reasonable approxima- 
tion, all relatives errors should be at least smaller than 
the unity, we assume that A < 1 from now on. As we 
have seen, the only non-trivial components of the Ein- 
stein equations that are exactly solved by the evolution 



law (5.3) are those corresponding to Gq and G\. For the 



remaining components, the error is at most of the order 
of I when the rotation parameter is small. 

Concerning our definition of relative errors, it is clear 
that e is the largest diagonal component of the energy- 
momentum tensor. In addition, we will see below that, in 
the spacetime region and range of parameters of physical 
interest, the other non- vanishing component of this ten- 



-Vo + V Vo^V 



V e [vi,vo] (■ ■ (6.1) 



In particular, for each fixed value of 77, the extrema of the 
interval Irf correspond to points that are causally con- 
nected with the observer. 

Moreover, the above region is invariant under the re- 
versal X — > —X. Using this fact and recalling that 
I € [0, 1) and D > 0, it is possible to show that, among 
all the conditions coming from the requirement that the 
relative errors be smaller than the quantity A, the most 
restrictive condition is that corresponding to the non- 
diagonal component Tq of the energy-momentum tensor. 
This component leads to the inequality 



le-'' GXia"^ + 6Da + 
Yi 3a2 -I- 6Da + 6A^ 



< A, 



(6.2) 



where a = a{ri) is given by Eq. ( ^.3| ), the pair of coordi- 
nates (77, x) must belong to the region ( |6.l| ), and we have 
adopted again the notation (5.4). 

On the other hand, from expression (|]^), we get 



2,0? + QDa + 6^2 ' 



(6.3) 



Recalling that the region under analysis is invariant un- 
der a flip of sign in the coordinate x (and that Xi > 0, 
Yi < 1 and A < 1), we then see that condition ( |6.2| ) en- 
sures that, in the region of physical interest, the energy 
density dominates over the non-diagonal component T2 
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of the energy-momentum tensor, as we had commented 
above. 

In addition, note that, since is a strictly decreasing 
function of x, its maximum value for a; G is obtained 
at (— ?7o + rf)/\/l — P . So, the most stringent condition 
contained in Eq. (|6.2|) is 



l<exp 



via) - yjao) 



3a2 + eDa + GA^ 



A. (6.4) 



We have employed here relation (5.6) to write the con- 
formal time in terms of a G [ai, oq], with oq > a\. These 
two values of the scale factor are reached, respectively, 
at the present time jyo and at the initial time of our con- 
siderations r]\. 

Using the explicit form of the function 77(a), it is actu- 
ally possible to show that, for fixed parameters / and A, 
the right-hand side of the above inequality is an increas- 
ing function of a. As a consequence, its minimum value 
in the interval [ai,ao] is attained when a — a\. In this 
way, we conclude that the necessary and sufficient con- 
dition for the relative errors to be smaller than A in the 



region of physical relevance is obtained from Eq. (6.4) by 
making a = a\. In particular, if we consider the whole re- 
gion that can be causally connected with the origin since 
the initial big bang, i.e. a\ — 0, we get 



l<exp 



3 A 

4^' 



(6.5) 



where we have employed that 77 vanishes when a = 0. 

This inequality sets an upper bound to beyond which 
our solution cannot be considered a good approximation 
modulo relative errors smaller than A. It is worth notic- 
ing that the conformal time ry(ao) that appears in Eq. 
( |6.5| ) depends on the rotation param eter as well as on 
the constants A and D, via relation (|5.6|). Owing to this 
dependence, it is in general difficult to find the exact 
value of the upper bound on / once the scale factor ag 
and the numbers A, and D are known. In the Ap- 
pendix, we present a method to estimate such an upper 
bound with great accuracy. In practice, nevertheless, it is 
possible to get a really good estimate by simply replacing 
Yi — \ — P with the unity and substituting "qia^) by the 
value 770 of the present conformal time corresponding to 
the exact solution with vanishing rotation parameter. It 
is not difficult to check that these approximations amount 
to disregarding corrections of the order of I- in the upper 
bound on I. Employing the values of ag. A, and D given 
in Sec. IV, one arrives in this way at 



I < 0.0337A. 



(6.6) 



As we have said, a more careful procedure to estimate 
this upper bound is presented in the Appendix, where we 
also consider the possibility oi = ao/1500, corresponding 
approximately to the time of decoupling between dust 
and radiation, and a model with slightly different cosmo- 
logical parameters, fid = 0.3 and H = 70 km/(sMpc). 



In all these cases, we obtain a value of the upper bound 
which is close to th e res ult given above. 

From inequality (3.6), we can easily derive an upper 
bound on the global angular velocity at present. Using 
that ui = l/{2a) and oq = 1.2x10^^ m (the value obtained 
in Sec. IV), we get w < 4.1 x 10~'^^A s^^ Thus, in 
order to have a relative error A < 2.5% one needs to 
impose, approximately, that u < 10~^^ s^^, while a more 
permissive error A < 25% would lead to a; < lO^^o g-i. 

Up to date, there exists no well-established and gen- 
erally accepted estimate of the angular velocity of the 
universe. In models with shear, some upper bounds can 
be inferred from the CMB anisotropy, but these bounds 
do not apply to the shear-free RTKO spacetimes. There 
are some estimations of uj based on the observed rotation 
of the plane of polarization of cosmic electromagnetic ra- 
diation leading to ~ 10^^* s^^. However, 
such observations are very controversial, and the derived 
value of Lo could well be two or three orders of magnitude 
smaller Q. 

An independent estimate uj ~ 10"^"'^ can be ob- 
tained from the analysis of the rotation of galaxies 
|30| . This result agrees with another estimation that is 
not based on observation, but on a heuristic argument, 
namely, the extension to the problem of rotation of the 
large number hypothesis put forward by Dirac. The an- 
gular momentum of the observed universe is L ~ puja^ , 
where p is the density of matter. From the large num- 
ber hypothesis, we get L ^ hA'jj where h is Planck 
constant and A^ ~ 10'^^ is Dirac scaling parameter js^. 
So, we have uj ~ hK%/{pa^). With p = 3 x 10^27 kg/m^ 
and a = ao = 1.2 X 10^^ m, this leads to w ~ lO"^! s"-"^. 

Let us finally remark that the upper bound that we 
have obtained for I is only aimed at determining the in- 
terval of rotation parameters in which the approximate 
RTKO solution presented in Sec. V is acceptably good. 
In principle, rotating solutions with larger angular veloci- 
ties are possible, but their energy-momentum tensor can- 
not be approximated by the matter content considered 
here. On the other hand, additional restrictions on the 
rotation parameter I could come from the requirement 
that the anisotropics that arise in the formulas of the lu- 
minosity and angular diameter distances are compatible 
with the observational data. The consideration of these 
anisotropies, however, cannot be carried out analytically 
if ? ^ 0, because, by contrast with the situation found in 
the case with vanishing rotation (see Sec. IV), the exact 
dependence of these distances with the redshift Z is not 
manageable anymore. What is available now is (the first 
terms of) their Kristian-Sachs expansion in powers of Z 
|25[ . Usin g th e expressions given by Obukhov for this 
expansion 124] and defining the relative variation of the 
angular diameter distance Sa like in Sec. IV, it is possi- 
ble to show that Sa — 21 up to second order corrections 
in Z and in the rotation parameter. Therefore, recalling 
the bound on I obtained above, we can affirm that the 
influence of rotation in the formulas for distances is negli- 
gible, at least as far as we do not consider sources of high 
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redshift. For high redshifts the Kristian-Sachs expansion 
is expected not to be vahd, and a more careful analysis 
is needed to determine the relevance of the anisotropics. 



VII. CONCLUSIONS 

In this paper, we have shown that it is possible to con- 
struct anisotropic models that are at least compatible 
with the main features of standard cosmology. In particu- 
lar, we have found an exact solution of Einstein equa- 
tions which describes an expanding universe containing 
an anisotropic scalar field and a comoving perfect fluid 
with two components: radiation and dust. The solution 
is spatially homogeneous, but the sections of constant 
time are anisotropic, its topology being the product of a 
pseudosphere and a real line. Even so, the background 
radiation is perfectly isotropic and the redshift experi- 
mented by any possible emission of the dust particles 
varies with the scale factor like in a FRW model. More- 
over, the expansion is conformal and follows the same 
evolution law as in a standard open FRW spacetime filled 
with dust and radiation. 

The relation between the redshift of astronomical 
sources and their angular diameter (or luminosity) dis- 
tance turns out to be anisotropic, because so is the space- 
time metric. However, this anisotropy does not conflict 
with the current observational data, because the corre- 
sponding variation of distances with the line of sight in 
our model is not dominant compared with the systematic 
and experimental errors of the measurements. 

The introduction of the massless, anisotropic scalar 
field leaves, nevertheless, one important imprint: the en- 
ergy density of the model is of the order of the critical 
density at all times. Therefore, the universe is always 
quasiflat. In more detail, the relative energy density 
equals the unity at the initial big-bang singularity, like 
in FRW cosmology, and decreases monotonically during 
the whole evolution to a lower bound of one third, which 
is the asymptotic limit reached at infinitely large times. 

The cosmological model that we have constructed is 
not completely realistic because, for instance, it does not 
predict the observed accelerated expansion of the uni- 
verse. In principle, this defect could be cured by in- 
cluding additional dark energy in the system, supplied 
either by a cosmological constant or by a quintessence 
field. This modification of our model will be discussed 
elsewhere. Here, we have concentrated our attention in 
our simple model because it permits a clear discussion 
of the effects of the anisotropic scalar field and allows to 
obtain explicitly the time dependence of the scale factor 
and the cosmological parameters. 

We have also presented a quite straightforward appli- 
cation of our exact solution, namely, the obtention of ap- 
proximate cosmological models describing spatially ho- 
mogeneous, anisotropic spacetimes with rotation. This 
has been possible because the anisotropic metric of our 



exact solution is in fact the element with vanishing ro- 
tation of a family of shear-free rotating metrics with re- 
markable properties, including the isotropy of the comov- 
ing CMB and the preservation of the standard relation 
between the redshift of light and the value of the scale 
factor when this light was emitted. 

Assuming that the matter content is the same as in 
our exact non-rotating solution, we have proved that it 
is possible to generalize the time dependence of the scale 
factor so as to attain an approximate solution of Einstein 
equations in the presence of rotation. More specifically, 
if one restricts all considerations to the causal past of 
the observer, we have shown that the error committed 
with our approximations in Einstein equations, relative 
to the energy density of the system (which is the dom- 
inant component of the energy- momentum tensor), re- 
mains smaller than any required quantity A if one sets 
an upper bound linear in A to the angular velocity of the 
present universe. In particular, we have calculated this 
bound using the values of the Hubble parameter and the 
relative energy density of pressureless matter provided 
by the concordance model |^^. For relative errors of a 
few percent, the upper bound that we have found turns 
out to be of the same order of magnitude as those ob- 
tained from observation of the rotation of galaxies | |30[ | 
and heuristic considerations involving the large number 
hypothesis 

Finally, an interesting possibility would be to analyze 
the angular power spectrum of primordial fluctuations in 
the CMB of these anisotropic cosmologies. This analysis 
would require an extension of the standard scheme of 
cosmological perturbation theory ||^ that dealt with the 
fact that the spatial sections of the spacetime are not 
maximally symmetric, took into account the anisotropic 
dependence of distances on the redshift, and treated the 
rotation parameter also in a perturbative manner. These 
issues will be the subject of future research. 
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APPENDIX 

In this appen dix, we will estimate the upper bound 
that inequality (|6.4|), evaluated at a = ai, sets to the 



rotation parameter I. Remembering expression (5.6), we 
can write the considered inequality as 



I 



{1-P) 



<J{Xi,Xi)A, 



(Al) 
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where 



J{U,V) = 



3al + 6Dai + GA^ 
6Uai + 6Dai + SA^ 



1 3Uai + D + ^W^a{ + WDax + 6^A^ \^^ 

Note that J depends on the non-negative constants A 
and D and on the values of the scale factor at present, 
ao, and at the initial time, a\. 

It is straightforward to see that, for oq > ai, J(J7, V) 
increases with V ^ assumi ng t hat U and V are positive. 
Since, according to Eq. ( |5.4|) , Xi ranges in (1/2, 1], it 
then follows that a necessary condition for Eq. (Al) to 
be satisfied is that //(I — Z^) < JiXi, 1)'^- In addition, 
one can check that IjiX — l^) is greater than J{Xi, 1) when 
I approaches the unity, whereas the opposite happens at 
I — 0, provided that ^ > 0. Therefore, the functions 
Z/(l — P) and J{Xi, 1) intersect each other at least once 
in / g [0, 1). Moreover, in this interval oil, both functions 
turn out to be strictly increasing. It is then possible to 
prove that the largest of the intersection points, 



L = max < / e [0, 1) 



I 



l-P 



J{Xul) 



(A3) 



can be obtained by numerical iteration. Namely, defining 
li = 1 and In+i = filn), one can get L as the limit of the 
sequence {In}, where f{l) = F[J{Xi, 1)] and 



F[J] 



VI + - 1 
2J 



(A4) 



Recalling then that A < 1 and 1 — P < 1, one easily 
concludes that a necessary condition for inequality (Al) 
to hold is 



I < J(Xl,1)A. 



(A5) 



L et u s now find a sufficient condition ensuring inequal- 
ity (|Al[ ). From our previous discussion, we already know 
that I < L and that J{U, V) increases with V ii U and 
V are positive. Employing the definition of Xi, we then 
see that J{Xi,Xi) > J{Xi,Xl). In addition, J{Xi,Xl) 
is an increasing function of I in the interval [0,1), re- 
gardless of the constant value of Xl G (1/2,1]. So, 
J{Xi,Xl) > J{1,Xl), since Xi becomes the unity at 
I — 0. Hence, it follows that a sufficient condition for 
Eq. (^) to hold is 1/(1 - P) < J{1,Xl)A or, equi- 
valently I < F[J{1, Xl)A]. Finally, taking into account 
that < J(1,Xl)A < 1 for aU the allowed values of A 
and Xl, and that F[J] > J(l - J^) if < J < 1, it is 
easy to derive the simpler sufficient condition 



I < J(1,Xl)[1- j2(l,Xi)]A. 



(A6) 



Using the values of oq, A and D obtained in Sec. IV, 
making ai = (i.e., considering the entire causal past of 



the origin) , and following the procedure explained above 
to determine the value of L, one can check that the ne- 



cessary and sufficient conditions given in Eqs. (A5) and 
(A6) lead in fact to coincident upper bounds on /, up to 
the third significant figure. With this degree of accuracy, 
one gets the bound I < 0.0337A, which reproduces in 
fact the estimate reached in Sec. VI. If one made instead 
ai = ao/1500, paying thus attention only to those events 
in the causal region which occurred (approximately) after 
the time of decoupling, one would obtain, with the same 
level of precision, / < 0.0442A. 

In order to check the sensibility of our estimates to 
the particular values adopted for the relative energy 
density of dust matter and the Hubble parameter, we 
have repeated the evaluation of the constants A and 
Z?, the present scale factor oq, and the upper bound 
on I taking fi^ = 0.3 and Hq — 70 km/(sMpc). In 
this case, following the arguments explained in Sec. IV, 
one gets A = 1.3 x lO^"* m, D = 7.2 x 10^^ m, and 
flQ = 1.1 X 10^^ m, which are close to the values found 
with Hq = 65 km/ (sMpc) a nd f id = 0.35. In addition, 
with fli = 0, Eqs. ( |A5| ) and (A6) lead now to the bound 
I < 0.0260A (again up to the third significant figure), 
whereas I < 0.0344A if oi = ao/1500. So, the upper 
bound reached for I is of the same order of magnitude in 
all the considered cases. 



G.F. Smoot et ai, Astrophys. J. Lett. 396, LI (19 92). 
For a recent review, see, e.g., E. Bertschinger, istro- 



ph/0101009 [Proceedings of COSMOLOGY 2000, edited 



by M.C. Bento, O. Bertolami, and L. Teodoro (in press)]. 
[3] P. de Bernardis et ai, Nature (London) 404, 955 (2000); 

S. Hanany et ai, Astrophys. J. Lett. 545, L5 (2000). 
[4] A. Balbi et ai, Astrophys. J. Lett. 545, LI (2000); A.E. 

Lange et ai, Phys. Rev. D 63, 042001 (2001). 
[5] R. Giovanelh et ai, Astrophys. J. Lett. 505, L91 (1998); 

D.A. Dale et al, ibid. 510, Lll (1999). 
[6] See, e.g., S. Weinberg, Gravitation and Cosmology (Wi- 
ley, New York, 1972). 
[7] D. Kramer et al, Exact Solutions of Einstein's Field 

Equations (Cambridge University Press, Cambridge, 

England, 1980). 
[8] V.A. Korotkii and Yu.N. Obukhov, Sov. Phys. JETP 72, 

11 (1991) [Zh. Eksp. Teor. Fiz. 99, 22 (1991)]. 
[9] J. Ehlers, P. Geren, and R.K. Sachs, J. Math. Phys. 9, 

1344 (1968). 

[10] A. Guth, Phys. Rev. D 23, 347 (1981); for a review see 
D.H. Lyth and A. Riotto, Phys. Rep. 314, 1 (1999). 

[11] C. Brans and R.H. Dicke, Phys. Rev. 124, 925 (1961). 

[12] A.G. Riess et ai, Astron. J. 116, 1009 (1998); S. PerL 
mutter et al, Astrophys. J. 517, 565 (1999). 

[13] J. Frieman et ai, Phys. Rev. Lett. 75, 2077 (1995); R.R. 
Caldwell, R. Dave, and P.J. Steinhardt, ibid. 80, 1582 
(1998); A.A. Starobinsky, JETP Lett. 68, 757 (1998) 



11 



[Pisma Zh. Eksp. Teor. Fiz. 68, 721 (1998)]; T.D. Saini 

et al, Phys. Rev. Lett. 85, 1162 (2000). 

L. Amendola, Phys. Rev. Lett. 86, 196 (2001); A. Balbi 

et al, Astrophys. J. Lett. 547, L89 (2001). 

J.D. Barrow and J. Magueijo, Phys. Lett. B 447, 246 

(1999). 

Yu.N. Obukhov, Gen. Relativ. Gravit. 24, 121 (1992). 
M.J. Rebougas and J. Tiomno, Phys. Rev. D 28, 1251 
(1983). 

W. Hasse and V. Perhck, J. Math. Phys. 29, 2064 (1988} 
S. Carneiro, Phys. Rev. D 61, 083506 (2000); 



qc/000309e 



gr- 



L.D. Landau and E.M. Lifshitz, The Classical Theory of 
Fields (Pergamon, Oxford, 1983), Sees. 111-114. 
S.W. Hawking and G.F.R. Ellis, The large scale structure 
of space-time (Cambridge University Press, Cambridge, 
England, 1973). 

W. Chen and Y.S. Wu, Phys. Rev. D 41, 695 (1990); 
M.V. John and K.B. Joseph, ibid. 61, 087304 (2000). 
See, e.g., M. Tegmark, M. Zaldarriaga, and A.J.S. Hamil- 
ton, Phys. Rev. D 63, 043007 (2001). 
Yu.N. Obukhov, in Colloquium on Cosmic Rotation, 
edited by M. Scherfner, T. Chrobok, and M. Shefaat 
(Wissenschaft und Technik Verlag, Berlin, 2000). 
J. Kristian and R.K. Sachs, Astrophys. J. 143, 379 
(1966). 

K. Tomita, Prog. Theor. Phys. 40, 264 (1968); M.A.H. 
MacCallum and G.F.R. Ellis, Commun. Math. Phys. 19, 
31 (1970). 

J.J. Cowan et al., Astrophys. J. 480, 246 (1997). 

R.G. Gratton et al, Astrophys. J. 491, 749 (1997); B. 

Chaboyer et al., ibid. 494, 96 (1998). 

R.W. Kiihne, Mod. Phys. Lett. A 12, 2473 (1997). 

L.X. Li, Gen. Relativ. Gravit. 30, 497 (1998). 

S. Carneiro, Found. Phys. Lett. 11, 95 (1998). 

P.A.M. Dirac, Nature (London) 139, 323 (1937). 



12 



